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ABSTRACT 


Statistical  analyses  of  the  dynamics  of  some  two-dimension¬ 
al  linear  vehicles  traveling  on  a  rough  track  are  performed  to  deter¬ 
mine  the  influence  on  two  aspects  of  vehicle  ride  of  a  set  of  parame¬ 
ters  which  include  wheel  base  length,  idealized  tire  imprint  length, 

•  peed,  and  damping  constant.  It  is  assumed  that  the  vehicles  move 
with  constant  horizontal  velocity  on  a  second  order,  weakly  station¬ 
ary  and  mean  square  continuous  random  track  with  contact  maintain¬ 
ed  at  all  times  between  the  idealized  tires  and  the  track.  The  two 
aspects  of  vehicle  ride  used  as  measures  of  the  ride  roughness  are 
peak  value  of  power  spectral  density  and  variance  of  frame  acceler¬ 
ation,  the  frame  acceleration  being  either  vertical  at  the  c.  g.  of 
frame,  vertical  at  the  point  over  idealized  wheel,  or  angular  (pitch¬ 
ing). 

For  the  same  speed,  damping,  power  spectral  density  for  the 
track,  and  two  particular  vehicles,  the  idealized  tire  imprint  length 
was  a  relatively  unimportant  parameter  over  a  fairly  large  range  of 
values.  On  the  other  hand,  one  parameter  which  included  the  wheel 
base  length  was  found  to  be  important  under  the  same  conditions. 

Four  sets  of  parameter  values  were  found  which  at  the  same 

•  peed  produced  best  or  optimal  rides  for  vertical  acceleration  at  the 
frame  c.  g.  and  over  the  wheel,  depending  upon  which  measure  of 
ride  roughness  was  employed.  The  influence  of  speed  was  then  ex¬ 
amined  on  vehicles  having  these  sets  of  parameter  values.  In  all 
cases,  increasing  speeds  produced  sharp  increases  in  ride  roughness 


ii 


PREFACE 

This  is  the  sixth  in  a  series  of  papers  studying  the 
motion  of  vehicles  under  random  excitation.  Prior  papers 
are: 

a.  Land  Locomotion  Report  No.  48,  Behavior  of  a 
Linear  One  Degree  of  Freedom  Vehicle  Moving  With 
Constant  Velocity  on  a  Stationary  Gaussian  Random 
Track . 

b.  Land  Locomotion  Report  No.  56,  On  the  Behavior  of 
a  Linear  Two  Degree  of  Freedom  Vehicle  Moving  With 
Constant  Velocity  on  a  Track  Whose  Contour  is  a 
Stationary  Random  Process. 

c.  Land  Locomotion  Report  No.  65,  On  the  Statistical 
Analysis  of  the  Motion  of  Some  Simple  Vehicles 
Moving  on  a  Random  Track,  which  also  includes  an 
appendix  which  defines  many  terms  and  derives  many 
basic  concepts  used  in  the  series. 

d.  Land  Locomotion  Report  No.  66,  On  the  Statistical 
Analysis  of  the  Motion  of  Some  Simple  Two-Dimensional 
Linear  Vehicles  Moving  on  a  Random  Track. 

e.  Land  Locomotion  Report  No.  72,  On  the  Statistical 
Properties  of  the  Ground  Contour  and  its  Relation 
to  the  Study  of  Land  Locomotion. 
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I.  Introduction 


The  motivation  for  the  investigation  reported  in  this  paper 
and  the  need  for  a  statistical  description  of  the  ground  in  the  inves¬ 
tigation  are  described  in  the  Introduction  to  the  preceding  paper. 

Hence,  the  remarks  of  this  section  only  need  be  directed  to  those 
aspects  of  vehicle  dynamics  and  "optimal  ride"  not  previously  com¬ 
mented  upon. 

A  wide  variety  of  one,  two,  and  three  dimensional  models 
have  been  suggested  for  a  study  of  the  dynamics  of  real  vehicles. 

Some  of  these  models  are  linear  while  others  are  non-linear.  The 
one  we  shall  employ  in  this  paper  [  1  ]  [  2  3  is  sufficiently  simple 
so  that  the  mathematics  is  reasonably  tractable  and  yet  sufficiently 
complicated  so  as  to  make  parameter  studies  interesting.  A  sketch 
of  this  linear  two-dimensional  vehicle  is  given  in  Figure  1. 

This  vehicle  consists  of  a  rigid  frame  connected  at  each  end 
via  a  parallel  arranged  spring  and  damper  suspension  element  to  an 

idealized  wheel-tire.  In  the  initial  stages  of  the  investigation,  a  num¬ 
ber  of  wheel-tire  models  were  considered.  However,  on  simplifying 
each  of  them  to  a  case  which  could  be  described  in  a  reasonable  math¬ 
ematical  form,  the  wheel-tire  model  shown  was  obtained.  It  has  mass 
but  no  moment  of  inertia,  and  it  is  connected  to  the  ground  or  track 
through  a  spring-damper  arrangement  distributed  along  a  massless 
rigid  bar.  Contact  with  the  track  over  the  entire  foot-pring  length  is 
assumed.  Thus,  the  idealized  tire  acts  much  like  a  deformable 
slipper  gliding  over  the  track.  Bouncing  is  not  permitted  in  the  sense 
that  contact  over  the  entire  foot- print  length  is  maintained  at  all  times, 
but  the  wheel  mass  may  oscillate  vertically  on  the  suspension  element 
and  the  idealized  tire.  We  also  assume  that  the  elements  shown  verti¬ 
cal  always  remain  in  that  position.  The  vehicle  of  Figure  1  thus  posses¬ 
ses  the  frame  translation  and  pitch  and  some  of  the  simpler  aspects  of 
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wheel-tire  dynamics  common  to  a  real  vehicle.  However,  it  is  clear 
that  frame  roll,  tire  bounding  off  the  track,  etc.  are  excluded. 

To  estimate  whether  one  vehicle  design  can  be  driven  faster 
than  another  design,  or  whether  a  given  vehicle  design  can  be  driv¬ 
en  at  a  definite  speed  over  a  terrain  of  specified  roughness  requires 
that  there  be  some  quantitative  ride  criterion  against  which  project¬ 
ed  vehicle  capability  can  be  measured.  At  the  present  time,  a  uni¬ 
versally  accepted  ride  criterion  is  not  available. 

It  is  a  commonly  held  view  among  individuals  associated  with 
the  military  wheeled-vehicle  field  that  such  vehicles  have  the  reserve 
power  and  sufficient  mechanical  strength  to  go  to  much  greater  speeds 
than  presently  being  attained  under  off-road  conditions.  Moreover,  if 
speeds  are  reached  that  begin  to  break  parts,  or  if  more  power  is  need¬ 
ed  to  increase  speed,  these  individuals  believe  that  the  needed  strength 
and  power  can  easily  be  supplied.  The  driver,  occupants,  or  cargo  of 
vehicles  thus  appears  to  provide  a  basic  limitation  on  vehicle  speed. 

The  frame  motion  of  the  vehicle  model  is  limited  to  vertical 
translation  and  pitch.  The  acceleration  of  each  appears  to  influence 
certain  aspects  of  ride.  For  example,  bouncing  of  the  driver,  say, 
upon  his  seat  is  associated  with  large  vertical  frame  acceleration. 

When  such  bouncing  becomes  sufficiently  severe,  it  results  in  physical 
injury.  Severe  oscillatory  acceleration  in  pitch  results  in  loss  of  vis¬ 
ual  perception  (this  also  occurs  with  severe  oscillatory  vertical  accel¬ 
eration)  and  loss  of  equilibrium.  In  most  cases,  the  driver  reduces 
vehicle  speed  before  either  the  oscillating  vertical  or  pitching  acceler¬ 
ation  reaches  such  extreme  values,  or  the  occupants  insist  that  he  do 
so.  Damage  tests  with  humans  to  determine  when  vertical  acceleration 

causes  injury,  pitch  acceleration  causes  loss  of  visual  perception,  etc.  , 
have  been  conducted  using  simple  harmonic  motions.  [3]  [4]  Results 
of  such  tests  cannot  be  directly  converted  to  the  case  in  which  the  oscil- 
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lations  are  random  in  character.  Hence,  bounds  on  vertical  accel¬ 
eration  (also  transverse,  etc.)  and  pitching  acceleration  set  by 
these  two  tests  are  not  useful  to  us.  Tests  in  progress  and  also  pro¬ 
jected  will  be  with  randomly  oscillating  acceleration,  and,  when  these 

results  are  available,  they  should  establish  quantitative  bounds  on  hu¬ 
man  tolerance  to  acceleration  in  random  oscillation.  In  any  event,  it 
is  reasonable  at  this  time  to  base  ride  criteria  on  random  oscillatory 
acceleration  in  vertical  acceleration  and  pitch. 

A  number  of  drivers  of  military  vehicles  also  pointed  out  that 
when  large  gun  tubes  are  carried  in  a  turret,  limiting  vehicle  speed 
was  frequently  set  by  violent  oscillations  occuring  in  these  tubes  as  a 
result  of  severe  y aide le  pitch. 

We  shall  select  two  quantities  as  measures  of  the  severity  of  a 
random  oscillation.  One  is  the  magnitude  of  the  largest  peak  in  the 
power  spectral  density.  The  smaller  the  maximum  peak,  the  lower 
the  ordinate  bounding  the  entire  spectral  density  curve.  While  the 
work  of  Goldman  C  3]  is  not  applicable,  it  does  suggest  that  low 
values  of  this  bounding  ordinate  may  be  associated  with  a  comfortable 
ride  from  the  human  standpoint  in  random  oscillatory  acceleration. 

The  variance  is  a  msasure  of  the  spread  of  a  probability  distri¬ 
bution  about  its  mean.  A  large  variance  implies  that  the  probability 
Of  a  significant  deviation  from  the  mean  is  large.  For  a  reasonably 
smooth  random  oscillation,  the  variance  may  thus  be  used  as  a  meas¬ 
ure  at  any  time  of  the  probability  of  the  oscillation  deviating  from  its 
mean  by  some  specified  amount.  Hence,  by  reducing  variance,  the 
probability  of  large  excursion  in  such  an  oscillation  is  in  general  re¬ 
duced. 

It  is  well  to  point  out  that  other  aspects  of  vehicle  acceleration 
may  provide  useful  ride  criteria.  In  particular,  the  work  of  von 
Gierke  [4l  and  others  [  5}  suggests  that  minimum  acceleration 
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power  spectral  density  {vertical  translation  or  pitch)  in  a  band  from 
3cpstol0cps  might  produce  the  best  ride  from  the  human  stand¬ 
point.  Thus,  the  peak  in  this  band  could  be  a  measure  of  the  sever¬ 
ity  of  a  random  acceleration. 

In  this  paper,  we  shall  study  the  influence  of  several  vehicle 
parameters  upon  an  optimal  ride  criterion  based  either  on  minimiz¬ 
ing  the  largest  peak  in  power  spectral  density  or  on  minimizing  the 
variance  in  vertical  translation  or  pitch  acceleration.  All  numeri¬ 
cal  results  presented  pertain  to  these  optimal  ride  criteria. 

We  shall  begin  the  main  portions  of  the  paper  by  presenting 
some  elementary  remarks  on  such  statistical  quantities  as  mean, 
variance,  covariance,  power  spectral  density,  m.  s.  continuity, 
etc.  as  they  apply  to  the  track  on  which  the  model  moves.  The  phys¬ 
ical  aspects  of  the  model  considered  in  the  equations  of  motion  are 
then  discussed.  Section  IV  presents  the  equations  of  motion  and 
shows  how  the  power  spectral  density  of  vertical  translation  of  the 
frame,  etc.  are  obtained.  Numerical  results  on  parameter  studies 
for  two  special  cases  of  the  vehicle  of  Figure  1  and  some  general  ob¬ 
servations  are  given  in  the  final  section. 


II.  Track  Model 


A  detailed  discueeion  of  a  stochastic  model  of  ground  rough¬ 
ness  is  being  presented  elsewhere  C  2D  .  It  suffices  for  the  purpos¬ 
es  of  this  paper  to  simply  assume  that  the  track  elevation  vs  horizon¬ 
tal  distance  is  a  random  function  possessing  statistical  properties  that 
ground  roughness  along  a  line  of  reasonable  length  is  likely  to  possess. 
In  particular,  we  shall  confine  our  assumptions  to  second  order  prop¬ 
erties  [7j  ,  i.  e.  to  properties  which  can  be  deduced  from  covari¬ 
ances.  Thus,  explicit  probability  distributions  6f  ground  roughness 
or  elevation  are  not  employed. 

The  mean  response  of  a  linear  vehicle  to  a  non- zero  mean 
track  elevation  may  be  treated  by  deterministic  and  well  known  meth¬ 
ods.  The  flucthations  of  the  response  about  the  mean  motion,  on  the 
other  hand,  must  be  treated  by  probabilistic  methods.  Since  we  are 
primarily  interested  in  the  latter,  we  shall  assume  that  the  mean  track 
elevation  has  the  constant  value  zero.  As  we  wish  to  employ  second 
order  properties,  it  is  required  that  the  track  of  elevation  be  a  second 
order  Q7  J  random  function  (  r.  f .  ).  Member  functions  of  ground  ele¬ 
vations  are,  of  course,  oscillatory  in  nature,  but  they  should  not  con¬ 
tain  discontinuities  involving  more  than  finite  jumps.  The  assumption 
of  mean  square  (m.  s.)  continuity  [7]  is  usually  sufficient  to  eliminate 
unreasonable  member  functions  and  we  shall  employ  it.  Finally,  in  or¬ 
der  to  be  able  to  employ  the  concept  of  power  spectral  density  and  the 
techniques  associated  with  it  C83,we  shall  assume  that  ground  elevation 
vs  horizontal  distance  is  at  least  weakly  stationary  [7J  .  This  means 
that  at  least  the  first  two  probability  distributions  do  not  depend  upon 
the  location  of  the  origin  of  the  horizontal  distance. 

As  an  additional  comment  upon  the  fact  that  results  given  in  this 
paper  only  depend  upon  second  order  properties  and  not  upon  particular 
types  of  distributions,  let  us  note  that  whether  a  particular  representa¬ 
tion  of  ground  elevation  as  a  r.  f ..  is  Gaussian  or  not  is  of  no  consequenc 
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There  are  essentially  two  ways  in  which  it  is  convenient  to 
represent  ground  or  track  elevation  Y0  (x)  as  a  r.  f.  of  the  hori¬ 
zontal  coordinate  x  .  In  one  way,  an  explicit  formula  for  Y0  (x) 
in  terms  of  a  family  or  set  of  random  variables  (r.  v's.  )  is  adopted, 
and,  in  the  other,  certain  statistical  requirements  are  postulated 
(zero  mean,  second  order,  m.  s.  continutiy,  and  at  least  weak  sta- 
tionarity)  and  a  certain  integral  representation  meeting  these  re¬ 
quirements  employed.  Each  procedure  has  certain  advantages;  hence 
we  shall  discuss  them. 

A  formula  adopted  in  previous  papers  [  1  ]  [  2  ]  is 


(2J) 


Y0  ( x )  *  oj  cos  (Aj  x  +  <pj) 


where  ay  and  Aj  are  ordinary  constants  and  the  tf>y  are  inde¬ 
pendent  r.  v's.  uniformly  distributed  over  an  interval  of  length  2  * 
Thus,  Y0  (x)  is  a  sum  of  cosine  functions  having  randomly  distributed 
phases.  The  Aj  have  the  dimensions  of  length  to  the  minuB  one  power; 
and,  when  divided  by  2  *■  give  cycles  per  unit  length.  Hence,  they 
may  be  regarded  as  angular  "length"  frequencies.  The  ay  are,  of 
course,  the  amplitudes  of  the  simple  harmonic  elements.  The  statis¬ 
tical  properties  of  (2.  1)  are  well  known  [  9  ] ;  we  easily  find  that 


(2.2) 


e  {ya  (*)}  =  o 

rYj*, ,  *,>  =  e  { w  s  <v}  -  2ji,  cos  aj  l*‘  ~  *'  ’ 


crj  (  x) 


n 


r  !x,*>  =  2 1 

Y0  jm! 


t, 


It  is  perhaps  appropriate  to  review  briefly  the  meanings  of  these 
equations.  The  first  states  that  the  mean  of  Y0  (x)  is  zero.  The 
third  states  that  fl-yf  (x)  is  independent  of  x  and  since  for  finite 
Os  and  finite  n  it  is  finite,  Y0  (x)  is  of  second  order .rY(x/txe) 
is  differentiable  any  number  of  times  with  respect  to  xf  and  xg  ; 
thus,  YQ  (x)  is  m.  s.  differentiable  any  number  of  times,  (In¬ 
deed,  the  member  functions  are  differentiable  to  any  order.  )  Final¬ 
ly,  since  fy  ( x,  ,  xg  )  depends  only  upon  xg-x,  and  the  mean  is 
constant,  we  know  that  Y0  (x)  is  at  least  weakly  stationary.  Hence, 
the  requirements  enumerated  above  for  the  random  function  which 
is  to  represent  ground  elevation  vs  horizontal  are  met. 

If  n  is  replaced  by  infinity,  there  are  conditions- -which  we 
need  not  state  here-  -  on  the  oj  and  \j  if  Y0  (x)  is  to  be  second 
order  and  m.  s.  continuous. 

The  power  spectral  density  of  Ya  ( x  )  is  given  by  the  formula 

[91  ' 

(2.3)  PYo(\>- 

where  8  (A  -  Xj)  is  the  delta  function  centered  at  Aj  with  the  proper 
ties 

f  S  (A  -  Aj)  dX  =*  i  ,  8  ( A  -Aj)  —  0  ,  if  A  i*  Aj 

-00 


We  shall  call  2  a*  the  weight  of  the  power  spectral  density  at  Aj . 

If  Y0  (x)  is  interpreted  as  a  voltage  applied  to  a  resistance  of  one 

ohm  and  x  is  interpreted  as  the  time,  then  the  integral  otPY(A)dA 

•o 

between  Aa  and  At,(>Aa )  is  the  average  power  in  Y0  (x)  between 
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those  values  of  A  .  On  setting  xg  -  x,  «  x0  we  may  write  down  the 
well  known  [  8  ]  formulas 


(2.4) 


Jy0  ~ 


y0 

t 


i 

I  Pv  (A) cos  A  xa  dA 

ro 

mOO 

PvJA)  -  jT  rYo  (xQ)  cos  A  xQ  dx0 

.00 

J  Py  ( AJdA 

Jn  rO 


°Y 
• O 


That  is,  the  covariance  /J,  ( xQ)  and  the  power  spectral  density 
Py  (x)  are  Fourier  transforms  of  one  another. 

If  a  set  of  the  <Pj  are  selected  with  the  aid  of  a 

table  of  random  numbers,  say,  then  one  member  function 


(2.5) 


n 


(x)  =  2 Z aj  cos  f  A/ 

/-/ 


*  *  #)> 


of  Y0  fxj  is  obtained.  Different  selections  of  the  j Aj  lead  to  different 
member  functions.  We  have  shown  elsewhere  [  1  ]  how  it  is  possible 
to  select  the  aj  and  Aj  so  that  the  integrated  power  spectral  density 

f  Pr  (A')dA’ 

J0  rO 

of  (2.  It)  approximates  as  closely  as  is  required  to  the  integral 

/ P( A')  d  A9 

Jo 

of  some  prescribed  or  given  power  spectral  density  P(A)  .  Since  sum 
a  procedure  determines  the  Oj  and  Aj  ,  it  is  thus  possible  with  (2.  ;>) 
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to  obtain  a  deterministic  member  function  y  (x)  of  the  random  func- 

o 

tion  (2.1),  and  the  integrated  power  spectral  density  of  this  YQ  (x)  ap¬ 
proximates  to  the  integral  of  the  prescribed  power  spectral  density 
P  (X)  .  A  graph,  or  for  that  matter,  an  actual  track*  of  a  section 
of  the  member  function  (2.5)  belonging  to  this  YQ  (x)  can  therefore  be 
constructed. 

The  advantages  in  employing  or  considering  (2.1)  are  that  a 
simple  and  readily  inter  pretable  formula  for  YQ(x)  is  available,  mem¬ 
ber  functions  are  easy  to  obtain  and  present  in  graphical  form,  if  desir¬ 
ed,  and  the  power  spectral  density  (2.3)  possesses  an  interpretation  in 
terms  of  well  known  electrical  engineering  concepts.  These  advantag¬ 
es  of  (2.1)  prompted  us  to  introduce  it  for  the  benefit  of  those  readers 
whose  background  in  random  function  theory  needs  refreshing.  The 
disadvantage  in  using  (2.1)  in  analytical  work  stems  from  the  fact  it 
often  leads  to  a  good  deal  of  unnecessary  algebra. 

Let  us  now  examine  the  formula 

(2.6)  Y0  (x)  =  faeiA*  </Z  (A) 

—  CO 

where  the  integral  is  to  be  interpreted  in  the  m.  s,  sense  [  7  ]  and 
the  random  set  function  Z  (X)  is  an  orthogonal  process  [7  ]  with 
the  properties 


A  track  calculated  on  this  basis  has  been  constructed  at  the  Detroit 
Arsenal. 


9 


(27)  1 


E  M  -  0 

jz£  {  1  *  <*>*)- i-  Py0<*> 


Py(X)=  Py  (-X) 
ro  To 


00 

T  J  py0  dA  <  oo  ,  -jf  J  / py  (A)  dA  < 

—oo  ° 


oo 


From  these  equations,  we  find  that 

r 

E  {yo(x)}-° 

,  OO 


(2.8)  ^ 


ry  ( x, ,  x£)  **  -j-  J  e /A  (  **-*/>  pYq  (A)  dA 


-  CD 
>  CD 


**  J  Py  (A)  cos  A  x0  dA 


CrY0(*)~~jrJ  Py  ( A  )  d A  <  OO 


-OO 
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Thus,  Yq  (x)  has  a  zero  mean,  is  of  second  order,  and,  since  the 
covariance  depends  only  upon  xg  —  x,  and  the  mean  is  constant,  it 
is  at  least  weakly  stationary.  The  last  of  (2.7)  and  second  of  (2,8) 
enable  us  to  demonstrate  the  m.  s.  continuity  of  (2.6).  Hence,  all 
the  requirements  which  we  wish  to  impose  on  the  r.  f.  representing 
ground  or  track  elevation  as  horizontal  distance  are  met. 

In  what  follows,  we  shall  employ  (2.6)  rather  than  (2.1),  be¬ 
cause  the  former  lends  itself  to  certain  computations  we  wish  to 
make . 
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III.  Vehicle  Model 


Figure  1  shows  the  idealized  linear  two-dimensional  vehicle 
that  is  to  be  discussed.  The  rigid  frame  has  wheel  base  length  l  ; 
its  mass  is  M;  and  its  moment  of  inertia  about  its  c.  g.  is  I  m .  The 
c.  g.  is  at  distance  b  from  the  wheel  base  center  O.  Y  ft)  is  the 
vertical  displacement  of  O,  and  9  ft)  is  the  angular  displacement  of 
the  frame.  We  shall  assume  that  Ohas  a  constant  horizontal  speed 
K  so  that  its  x-coordinate  is  vt,  t  being  the  time. 

The  ends  of  the  frame  are  supported  on  linear  springs  and 
viscous  dampers  arranged  in  parallel.  kj  and  Cj  (  j-!fl-) are  the 
rear  (j-U  and  front  (  j-2)  spring  and  viscous  damper  constants  re¬ 
spectively.  These  elements  are  taken  as  massless,  the  upper  and 
lower  pivots  are  assumed  frictionless,  and  they  remain  vertical 
throughout  the  motion. 

Idealized  tires  are  assumed.  Each  has  a  rigid  massless 
frame  of  length  2 aj  and  a  mass  mj  concentrated  at  the  frame  cen¬ 
ter.  These  frames  are  connected  by  frictionless  pivots  to  the  bot¬ 
toms  of  the  parallel  arranged  spring  and  damper  elements,  as  shown. 
Kjf(j)  and  rjffj)  are,  respectively,  the  linear  spring  constant  per 
unit  length  and  the  viscous  damping  constant  per  unit  length  which  are 
distributed  in  a  parallel  arrangement  over  the  bottoms  of  the  frames; 
we  assume  these  elements  are  massless.  We  shall  make  the  initial 
assumption  that  Kj((j)  and  ^  (£])  are  even  functions  of  (j  ;  at  a 
latter  stage,  we  shall  simply  assume  they  are  independent  of  (j  . 
Contact  of  the  tire  spring  and  damper  elements  with  the  track  a  road 
is  assumed  to  be  maintained  at  all  times.  9,(t)  is  the  angular  dis¬ 
placement  of  the  rear  tire  frame  and  9g(t)  is  the  angular  dispplace- 
ment  of  the  front  frame.  The  vertical  displacements  of  the  masses 
mj  and  m  2  are  Yt(t)  and  Yg  ft),  respectively. 

Gravity  will  be  neglected.  Further,  all  angles  are  assumed 
to  be  sufficiently  small  so  that  cosines  may  be  replaced  by  unity  and 
the  sines  by  the  angles  in  radians. 
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IV.  Analysis  of  Vehicle  Dynamics 


For  the  vehicle  model  just  described,  we  find  the  equations 
of  motion  to  be 


M  (Y  +  b$)  =  -F'(t)-  Fg(t)  , 


(4.0 


(-j-  +  b)F,(t)-  (-j--b)  Ft  (t)f 
m,  Y,  -  Ft(  t )  —  Fs(t)  t 


mt  Yg  -  Fg(t)-F4(t) , 


where 


Fl(t)-k,(Y-±®-Yl/+cl(Y-j-®-Yl}  , 
Ft(t!-k,(Y+±@-  Y,)+ct(9+±  &-Y,)  , 


14.2/  i 


°t 

+  !  <f,> [*,+(,  &, ~ *. 

F4(t)=f°‘«,(('){Y'+('lSl'-Y0(,+±+Si] 

~°g 

+z<fAi+$4*-i'o(*+i+fA}d<, 


and  where  derivatives  and  integrals  are  to  be  interpreted  in  the  tn.  s. 
sense  [7  ] . 
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To  simplify  the  analysis  from  this  point  on,  it  is  convemeni 
to  assume  that  the  vehicle  is  symmetrical,  i.  e. 


(4.3) 


k,  =  ke  =  k  ,  c,  —  ct—  c  , 

o,  =  og=o  t  (£ )  =  ^g( € )  =  k(£)  » 

r,  (£)=%(£)=  X(() ,  m,  —  mt  =  m  , 

b  =  0 


The  equations  of  motion  may  now  be  written  in  the  form 


r  .. 


Y+2£y  OJy  Y+  OJy  Y-  -£r(2Cy  dJy  Y'+  UJyYl)-j(2£y  (JyYg+vf  Yg)  -O 

§+2[  aJe®+6L>g0+i(2La>ei+vJr,>-{(2Ct“A+a'»  Yt>  =  0 


(4.4)  < 


- Y  +  ojSy)*^-  (&g*-®+v‘o®) 

+  Y,  +  (2(  Ofc*  <t,  —&~) Yt  +  a,o  <!+«,>*, 

*oW, 

y+uiy)-  (*%%*  -&+«>;  ®) 

+  Yt+ (zCouo+ai  *^°~-)y*  +  wZ  <l+aCi>Y‘ 

=  U,; J° [  Y0(*+±+(>+  -*£-  *.(>*  +  *f) j 
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where  now 


(4.5) 


,  2k  o7  -  2c_  ... /  k  l* 

-  ~ar  •  m  *  °b-  -jj—  • 

zifaJe=‘~§i^  •  a,‘~~w  •  zC0wo^  "w  • 

a,  -  -f  .  *-!'*  <(>*(  .  c  ~J°r(Odf 


These  symbols  have  the  usual  interpretations  given  in  small  vibra¬ 
tion  theory.  The  angles  0t  ( t)  and  0g  ( t)  drop  out  of  the  third 
and  fourth  of  (4.  4)  since  we  have  taken  K  ( ()  and  r(()  to  be  even 
functions  of  /• 

We  must  now  employ  a  representation  for  Y0(x) .  As  mention¬ 
ed  above,  it  is  convenient  to  use  (2.  6) .  It  is  also  convenient  at  this 
time  to  assume  the  forms  of  K  (()  and  Y( ()  are  the  same.  Thus, 

M(l  m  wfj.  m  lill  , 

K  C  2a  ' 

where  g  ( ()  is  even  in  (  and  the  two  non  zero  terms  on  the  right 
hand  sides  of  (4.  4)  become 

(4.6)  G(^)(l+'%i)  eiwt  dz(y-) 

-00 

respectively,  where 


(4.7)  cu —A  v  ,  G  (-~)  -  -jj  f  g(()eU¥~  d( 

-a 
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With  v  a  constant,  YQ  (vt)  and  Y0  (vt)  are  still  at  least 
weakl .  stationary.  It  follows  [  8  ]  that  Y(t)  ,  B-(t)  ,  Yt(t)  ,  Yg(t) 
are  also  at  least  weakly  stationary,  since  they  are  related  to  YQ  (vt) 
and  its  m.  s.  derivative  Y0  (vt)  through  a  system  of  constant  co¬ 
efficient  linear  differential  equations.  Thus,  we  ma\  write 

t  00  iCUt  ...  ,  , 

J  e  dU(cu) 

-00 

[°°  e,(ut  dV(o))  ' 

-CO 

f°°e  i0J1  dW,(w) 

-CD 

Yt(t)  =  J°°e  /aft dWg(oj) 

_  -00 

the  random  set  functions  U  (*>)  ,  V  (•»),  W,  (•*),  %(»)  being  orthogonal 
processes.  The  substitution  of  (4.  6)  and  (4.  d)  into  (4.  4)  provides 
us  with  the  following  equations  for  the  determination  of  dU  (*)  , 
dV(*>)  ,  dW/M  and  dWg(<*)  1 


(4.8) 


Y  (t)  * 
®(t)  = 
Y,(t)  = 
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(4.9)  < 


(cu*-  tu*+  i2 £  cuyCu)dU(cu) -  -j-  (cu*-hi2^  cuy  cu)dWt  (cu) 
-  4r  (ojy  +  i2  £  cuy  cu)  dWg  (cu)  =  O  , 

(ajg  -  ou‘+  i2£g  ojgQj)  dV(cu)+  (oj*  +  i  2 £#  cu&  cu)  dW,  (cu) 

--L-  (aj*+i2£ffcugcu)dwt(cu)  =0  , 


-  ac,  (oj*+  i  .*£x£jB—)dU(cu)+  -^1  (aj‘+i  )dV(cu) 

UJy  U/g 

+  [  OJ*  d+ac,)-  ou2+  i( ?£o  cu0cu  +  ac,  .2£y“<?OJJ  ]dW/cu) 

U/y 
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Tlie  solution  of  the  first  two  of  (4.  9)  for  dW,  (uj)  and 


gives 

r 

d  Wt  (cv)  =  — 

w,dUfco)  +  WgdV  (co) 

(4J0) 

w 

dWa  (cu)  =r  - 

_wj  dU(co)  +  w.dVfQj) 

- 

w 

with 

(4./1) 


"  =  (Vy*  i?£y  coycuj  (co£+i2^  c^cu)  , 

w<  ®  2  (uj0  +> 2 Vguj)( at*  -  ou‘+  i2£y  ojy  a>) 

w*~  l  (ojy-b i2{  ujy cu)( cuj-cu*  +  i2£  oj0cu) 
y 

The  last  two  of  (4.  9)  with  (4.  10)  determine  dU(cu)  and  dVfco): 
dU(cu)  — 

dV(cu)  = 


(4J2) 


where 


(4.13) 


- - cos  -§^-  dZ  (—)  , 

a.  +  Jlgi.  gv  x  v  '  ' 


#,  ,  w 

—iz 

v  +  w 


sin  dZ  (-%-)  , 


Ul=ai(U'+il£^L)  , 

$ 

u  ~  tai‘  a+cc,)~  We+i  ( ?£o  cu0cu  + a, 

OJy 

2  =  a,‘  G(^S-)  (l  +  i 

r  1  cun  / 


$>  fou; 
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The  power  spectral  density  of  Y(t)  ,  for  example,  is  given 
by  a  formula  similar  to  the  second  of  (2.  7);  it  is 

(4.14)  -TPy 

Hence,  with  the  first  of  <4.  12)  and  the  second  of  (2.  7),  we  find  that 


(4.15) 


*s  *'*  * a> 

+  /J&SLJf]"(i+cos  *f)  [e 


and  similarly. 


-  ("«-  $  -4f ^  ('-■$ 


2 


+  / 


^  I ')]"«  -cos -If)  [fi(-e*)]\  (W 


Two  special  cases  of  the  system  shown  in  Figure  1  are  of 
particular  interest  in  this  paper.  The  first  is  obtained  by  replacing 
the  idealized  tires  by  massless  point  follower.  The  power  spectral 
densities  of  the  vertical  displacement  Y  (t)+  z  B(t )  at  a  point  on 
the  frame  at  a  distance  z  along  the  frame  from  the  c.  g.  and  of 
B(t)  are 


The  power  spectral  densities  of  Y,(tJ,  Yg( t )  and  of  the  vertical  motion 
at  points  on  the  vehicle  frame  other  than  its  c.  g.  may  be  found  in  the 
same  manner. 
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respectively,  where  the  frequency  range  is  from  —  00  to  00  • 

The  second  case  is  obtained  from  the  system  shown  in  Figure 
by  making  the  two  elements  connecting  the  frame  and  the  tire  mass¬ 
es  rigid  and  lumping  the  tire  masses  in  with  the  mass  of  the  frame; 
we  find  that  the  spectral  densities  of  Y  (t)  +  Z  0(t)  and  0  (t) 
are 


(4.18) 


Py(z,(t))t  **  [  G  P S  (z  l  &),  , 

r 

Pe(z.u>),=  ■ 


We  have  omitted  derivations  of  (4.  17)  and  (4.  18)  since  they  follow 
almost  the  same  sequence  of  steps  that  we  employed  in  obtaining 
(4. 15)  and  (4.  16)  from  (4.  1 ). 
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V.  Discussion  and  Results_ 

It  is  possible  to  draw  a  few  general  conclusions  from  the 
results  given  in  the  previous  section. 

Examination  of  (4.  15)  and  (4.  16),  say,  reveals  that  the  half¬ 
tire  frame  a  half  foot  print  length  a  only  enters  these  equations 
through  the  term  G  (St^L)  where  [  see  (4.  7)  ], 

-a 

and  9(()  is  the  distribution  along  the  frame  of  the  parallel  ar¬ 
ranged  linear  spring  and  viscous  damper  tire  elements.  Since  6(0?-) 

iajf  v 

is  an  integral  of  g  (()  #  \t  /2a  over  the  foot  print  length, 

G  ( -^p~)  represents  the  smoothing  of  the  roughness  in  the  track  due 
to  the  idealized  tires.  Thus, 

(5.n  [e 

may  be  regarded  as  the  power  spectral  density  of  the  input  to  the  ve¬ 
hicle  obtained  from  the  power  spectral  density  of  Y0(x)  after  smooth¬ 
ing  by  the  idealized  tires.  This  suggests  that  in  the  analysis  we  could 
have  ignored  the  tire- foot  print  length  provided  the  track  Y0  ( x)  with 
power  spectral  density  Py0(~^v~)  were  replaced  by  a  smoothed  track 
with  power  spectral  density  (5.  1).  Advantage  of  this  point  may  there¬ 
fore  be  taken  in  future  investigations.  The  conclusion  that  a  tire  is 
a  smoothing  device  in  some  sense  is  in  agreement  with  re  suits  ob¬ 
tained  on  large  low  pressure  tires. 
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In  the  numerical  results  presented  below,  we  shall  assume 

t  hat 

(5.2)  g  ({)  =  /  , 

and,  thus, 

r  .sin  0 

(5.3)  [  G  (-&?-)]  =•(■  ya  ) 


It  is  clear  from  this  equation  that  the  larger  a,  the  narrower  the 
peak  in  the  term  on  the  right.  Hence,  the  larger  a,  the  smaller 
the  high  frequency  portion  of  (5.1)  and  the  greater  the  smoothing 
influence  of  the  tires.  Conversely,  the  smaller  a,  the  broader  the 
peak  in  the  right  of  (5.3),  and  the  less  the  smoothing  effect  of  the 
tires.  Actually,  the  selection  (5.2)  is  not  particularly  realistic, 
since  g(()  should  approach  zero  as  ±a.  However,  (5.3)  is  easy 
to  compute  with,  provides  an  estimate  of  the  influence  of  the  tire  on 
spectral  density  of  various  quantities,  and  so  is  adequate  for  pres¬ 
ent  purposes . 

The  power  spectral  density  Py(coJ  con.ta.ins  the  wheel  base 

length  l  in  just  one  term,  namely  in  (I  4-COS  COl/v  ASince  this  term 

multiplies  Py  (cj/v)  ,  just  as  does  [6  ( co  a/v  )~\ ,  it  may  also  be  re- 
ro 

garded  as  a  filter  acting  on  Py  ( <X)/v ) .  Clearly,  the  term  filters 
periodically,  depending  upon  the  value  of  l/v  .  Thus,  wheel  base 


length  and  speed  v  may  be  expected  to  have  a  substantial  influence 
on  certain  aspects  of  vehicle  ride.  The  te rm  ( l-COS  OJlAf)  in  P^(co) 
may  be  interpreted  in  the  same  sense. 
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From  the  first  of  (4.8),  we  deduce  that 


(5-4j  rr(, ,,,,!=!■  f°e  Py  lw)dw  . 

—CD 

The  r.f.  Y(t)  has  two  m.  s.  derivates  t  (tJ  and  Y(t)  •  it  follows 
[8]  that  the  covariances  of  these  r.  f.  are 


(5.5) 


* 


Hence,  on  combining  (5.4)  and  (5.5)  and  making  use  of  the  definition 
of  the  covariances  of  f  (t)  and  Y  (t)  in  terms  of  integrals  of  their 
power  spectral  densities,  we  find  that 


(5.6) 


Pf(Oj)  -  OJ*  Py(Cd)  , 

* 

Py(aj)~  cu4Py(cu)  ; 

<_ 


thus,  the  power  spectral  density  of  the  first  or  second  m.  s.  deriva¬ 
tive  of  a  suitable  r.  f.  may  be  found  by  multiplying  its  power  spectral 
density  by  or  by  cu4  ,  respectively.  This  rule  is  quite  general  C 8 J 
and  we  shall  use  it  without  further  comment  in  what  follows  to  obtain 

•I  „  •• 

the  power  spectral  densities  of  9  (t)  ,  and  Y  (t )  +  z  9  (t) 

The  constants  which  define  the  vehicle  and  track  enter  the  spec¬ 
tral  densities  of  interest  in  a  fairly  complex  manner.  Further,  the 

measures  of  ride  roughness  that  are  being  employed  involve  these  pow¬ 
er  spectral  densities  either  through  an  integral  (variance)  or  a  deriva- 
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tive  (peak  value).  Hence,  we  shall  not  employ  analytical  methods 
to  perform  quantitative  parameter  studies  on  ride  roughness  but  in¬ 
stead  use  numerical  methods . 

A  definite  formula  for  PyfXj is  now  required 

'O 


(5.7) 


Py  (A)  - 


AoVZir 


The  quantity  is  the  variance  of  YQ(x)  /  \Q  is  a  constant  which 

may  be  used  to  adjust  the  spread  of  the  single  peak  possessed  by 
PYo(X).  This  formula  appears  to  give  a  reasonable  fit  to  a  power 
spectral  density  which  we  obtained  from  ground  measurement.  Fur¬ 
ther  measurements  may  well  suggest  different  expressions  for  the 
right  hand  side  of  (5.7)*.  Hdwever,  for  the  present,  it  is  a  reason¬ 
able  selection. 

A  set  of  dimensionless  parameters  we  have  found  convenient  to 
employ  in  numerical  calculations  is  the  following: 


(5.8) 


a  «r 


T-  ’  4 


where  vQ  is  a  reference  speed  and  p  is  the  radius  of  gyration  about  r 
the  c.  g.  of  the  frame.  In  terms  of  these  parameters,  the  first  of 
(4.17)  can  be  rewritten  as 

*  The  covariance  corresponding  to  (5.7)  i6  analytic;  thus,  the  r.  f. 
having  this  covariance  is  m.  s.  analytic  [8  Ch.  7]]  While  m.  %. 
analyticity  has  certain  predictibility  implications,  we  shall  ignore 
them. 
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<ss)  V'-#r/]  {" 

y 

+  cosCC  l/r)\l+ ( 2£  r)*\\(l- ~=t)  +  (2  £yr)  ] 


*  8Cy  '*  T  sin  *  l/r  \(w  ~  /  “  4^r‘)~  ~^T (r-l-4C‘yr*)  ] 

+  *7T  (l-cosecrr)[t+(2^  ]} 6 


for  example.  The  other  spectral  densities  can  be  put  in  a  similar 
form. 

As  mentioned  in  the  Introduction,  the  two  ride  criteria  which 
we  shall  employ  are  either  the  magnitude  of  the  largest  peak  in  the 
power  spectral  density  or  the  value  of  the  variance  for  acceleration. 

t#  .  . 

The  accelerations  of  the  frame  we  shall  consider  are  Y  (1)  , 

Y  (t)  +  -£  •  (t)  and  O  (1)i  the  first  is  vertical  acceleration  at  the 
c  .  g.,  the  second  is  vertical  acceleration  over  idealized  wheel-tire, 


and  the  third  is  pitching  acceleration. 

Only  the  two  special  cases  of  the  vehicle  shown  in  Figure  1 
will  be  considered  at  this  time.  The  formulas  needed  for  the  power 
spectral  densities  are  obtained  from  either  (4,17)or  (4.18).  Ordin¬ 


ate  values  for  vertical  acceleration  results  are  divided  by  crJ  C*)^ 

'0  * 

and  those  for  angular  acceleration  results  by  Oy  these  factors 

will  be  omitted  on  the  graphs.  We  shall  also  omit  the  subscripts  1 


or  2  on  the  ordinate  designation,  sinc  e  the-  subsc  ripts  I  and  2  re¬ 


fer  to  the  cases  a/l  —  O  and  o/l  ^  0  ,  ri  spectively,  and  these  are 

indicated  explicitly  upon  the  graphs. 
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Figures  2-4  present  numerical  results  on  the  maximum  value 
of  the  power  spectral  densities  and  on  the  variances  of 
Y(t)  ,  Y(t)+Z0(t)  and#  (t)  as  functions  of  the  parameters  <X  and  o/l 

It  is  immediately  clear  that  results  on  variance  with 
o/l  *0  and  a/l  *  0.0625  and  also  on  magnitude  of  the  maximum  val¬ 
ue  of  the  power  spectral  density  do  not  differ  appreciably  until  ac  be¬ 
comes  larger  than  4  or  5.  We  note  that  o/l  ■  0.0685  cor  res  ponds  to 
a  total  foot  print  length- -rear  plus  front--  of  25%  of  the  wheel  base 
length  If  most  wheeled  vehicles  have  a  value  of  approximately  10%. 
Thus,  a/l  does  not  appear  to  be  a  significant  parameter  with  the  ride 
criteria  we  are  employing  unless  ac  exceeds  4  or  5.  In  what  follows, 
therefore,  it  is  convenient  to  set  a/l"*0  . 

On  the  other  hand,  changes  in  the  parameter  ac  have  a  pro¬ 
nounced  influence  on  the  magnitude  of  the  largest  peak  in  the  power 
spectral  density  and  the  variance  of  the  various  accelerations.  This 
strong  influence  is  due  to  the  fact  that  the  location  of  the  bands  on  the 
co  -axis  which  the  trigonometric  terms  containing  OC  filter  change  with 
cc  ,  sometimes  reducing  peaks  sharply  and  sometimes  reinforcing 
them. 

For  the  vertical  acceleration,  there  is  a  unique  value  of  ac  which 
makes  the  peak  magnitude  and  variance  a  minimum.  Thus,  if  the  best 
ride  under  the  conditions  assumed  is  determined  either  by  the  minimum 
peak  magnitude  in  the  spectral  density  or  by  the  minimum  value  of  vari¬ 
ance,  we  see  that  for  each  point  on  the  frame  there  is  a  unique  value 
OC  which  provides  such  a  best  ride;  the  range  of  cc  values  being  from 
.4088  to  2.6150.  We  also  note,  however,  that  slight  changes  from 
these  CC  values  produced  large  changes  in  minimum  peak  magnitude  and 
variance,  i.  e.  the  best  ride  is  a  sharply  tuned  function  of  cc  .  For  too 
large  a  value  of  OC  with  Z  =*  the  variance  is  larger  than  with  a  very 
small  OC  ;  but  with  Z  =  O  ,  the  variance  is  always  reduced  below  its 
value-  at  oc—O  by  increasing  OC 
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The  peak  value  in  the  spectral  density  and  the  variance  for  the 
pitching  acceleration  exhibit  a  difference  dependence  upon  a  than  en¬ 
countered  with  the  vertical  acceleration,  as  may  be  seen  from  Figure 
4.  Each  is  zero  when  OC  equals  zero;  as  ac  increases  each  reaches 
a  maximum  value,  falls  to  a  minimum  and  then  oscillates.  The  value 
of  OC  which  gives  the  best  ride  in  vertical  acceleration  does  not  give 
a  best  ride  in  pitch.  However,  it  must  be  remembered  that  the  param¬ 
eter  S  is  adjustable.  Hence,  it  may  be  possible  to  select  oc  and 
g  to  give  a  best  ride  from  the  point  of  view  of  vertical  acceleration 
and  pitch,  or,  at  least,  a  good  compromise  to  a  best  ride  in  both 
cases;  this  point  is  now  under  investigation. 

jf.y  *¥  • 

The  ride  criterion  used  in  Figure  5  is  the  magnitude  of  the  larg¬ 
est  peak  in  the  power  spectral  density  of  9  ( t )  .  With  —  *  0  »  fl—  1.00 
it  am  1.00  f  the  value  of  OC  which  makes  the  peak  a  minimum  is  plotted 
as  a  function  of  Cy  (See  Figure  5a).  In  addition,  the  minimum  peak 
value  is  plotted.  Although  the  range  of  the  minimum  peak  is  not 
large  for  the  considered,  there  is  a  unique  Cy  (»  .6953)  and 
oc  («  1.9170)  which  makes  this  minimum  peak  a  minimum.  Thus, 
one  pair  of  Cy  and  OC  give  the  best  ride  in  the  peak  of  power  spectral 
density  sense.  Figure  5b  is  a  graph  of  the  peak  in  the  power  spectral 
density  of  V(t)  as  a  function  of  r  with  -j-  =  O  ,  1.00,  ^*.6953, 

oc  ■»  1.9170.  A  r>/  corresponds  to  v<v0  and  *  <  /  corresponds  to  ¥>v0. 
Thus,  as  speeds  decrease  below  v0  ,  which  may  be  looked  upon  as 
design  speed,  the  peak  value  falls  gradually;  however,  as  k  increas¬ 
es  above  vQ  the  peak  value  rises  sharply,  i.  e.  the  ride  rapidly  gets 
worse . 

The  results  presented  in  Figure  6  correspond  to  results  given 
in  Figure  5  except  that  the  ride  criterion  is  now  the  value  of  the  vari¬ 
ance  of  9(1).  Figures  7  and  8  correspond  to  Figures  5  and  6  but  for 
9  (t)+  -g-  &  (t)  .  The  conclusions  drawn  from  Figure  5  are  reflected 
in  these  three  Figures. 
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Figures  8  and  9  present  resuits  on  0  (t)  that  is  the  saint 
as  given  in  Fig  ures  5  and  6,  respectively,  for  Y  (t)  .  However, 
with  0(t)  neither  the  minimum  variance  nor  minimum  peak  value 
has  a  minimum  value  as  a  function  of  Cy  -  This  result  agrees 
with  the  observation  that  the  oscillations  of  6  (t)  may  be  reduced 

to  zero  by  making  l  infinite  and  using  a  rigid  suspension  system. 

For  arbitrarily  selected  Cy  ,  these  two  Figures  again  show  the 
great  sensitivity  of  variance  and  maximum  peak  value  of  power 
spectral  density  to  changes  in  r  in  the  neighborhood  of  *»  1 . 

The  parameter  fi  was  not  changed  in  the  studies  reported 
in  this  paper.  We  recall  that  fi  contains  XQ  which  is  a  measure  of  the  s  pread 
of  PYJx)  .  At  this  time,  we  have  no  numerical  value  for  X0 
moreover,  we  are  not  certain  that  (5.7)  is  the  most  reasonable 
choice  for  PYo  (X)  .  Hence,  it  does  not  seem  appropriate  to  obtain 
results  with  other  than  one  value  of  fi  .  It  is  worthy  of  note,  how¬ 
ever,  that  since 

(5 JO)  CC/&-  X0  1  • 

if  we  had  a  numerical  value  for  X0  and  the  ac  associated  with  a  best 
ride  in  same  sense  for  a  given  Cy  t  &  *  etc. ,  then  a  numerical  val¬ 
ue  for  the  corresponding  wheel  base  length  could  be  obtained. 

The  results  presented  in  this  paper  represent  the  beginnings 
of  an  attack  on  the  interesting  problem  of  how  to  design  a  vehicle  to 
achieve  higher  speeds  than  obtained  at  present  under  off-road  condi¬ 
tions.  A  large  effort  will  surely  have  to  be  devoted  to  obtaining  a 
statistical  characterization  of  ground  contour  geometry  and  quantitative 
ride  criteria  before  either  the  method  of  attack  proposed  can  be  justi¬ 
fied  here  and  definitive  design  results  can  be  obtained. 
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